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Abstract 

We study the heat kernel for an operator of Laplace type with a 5- 
function potential concentrated on a closed surface. We derive the general 
form of the small t asymptotics and calculate explicitly several first heat 
kernel coefficients. 



1 Introduction 

t> , 

Singular potentials are a frequently used idealization of physical situations allow- 
ing for an easier (and, sometimes, explicit) solution while keeping the essential 
features of the problem. The best studied cases of singular potentials are the 
delta function potentials concentrated at isolated points, which describe contact 
interactions of particles (for a review, see |]). Rigorous analysis of such potentials 
was initiated by the paper by Berezin and Faddeev || and has developed later in 
<D \ a mature mathematical discipline ||TJ. Other cases of singular potentials include 

also cosmic strings and other topological defects || and problems related to the 
black hole entropy H. We like also to mention a recent work on the boundary 
discontinuities [§[]. 

With respect to the Casimir effect, a delta function shaped potential provides 
the simplest generalization of the conductor boundary conditions towards inclu- 
sion of more realistic properties of the walls like partial transparency. For a scalar 
and a spinor field with plane boundaries this problem has been investigated in 
@, for moving partly transmitting mirrors in 0. An interesting approach us- 
ing 'semihard' and 'weak' boundaries is developed in |J. In all those cases it 
is crucial to know the ultra violet divergences in order to find out the structure 
of necessary counterterms. This is equivalent to the investigation of the corre- 
sponding heat kernel asymptotics. As it is known ||, this is an expansion with 
respect to integer powers of the proper time parameter t (see below) for Laplace 
type operators on closed manifolds and to half integer powers, i.e., to powers of 
y/t, on manifolds with boundaries and Dirichlet or Robin boundary conditions. 
For more complicated pseudo differential operators powers of hit may appear. 
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The present paper is devoted to singularities which are located on closed hy- 
persurfaces of dimension m — 1 where m is dimension of the underlying manifold. 
Apart from the quantum mechanical problem of a particle in the space with 
semitransparent boundaries (for a recent calculation of the vacuum energy for 
such system see e.g. ||10| ) possible physical application also include fermions on a 
background of a magnetic tube [ljj and photons interacting with dielectric bodies 
U, HI- In all these cases the dynamics is described by a second order differ- 



ential operator of Laplace type supplemented by certain matching conditions on 
a surface. Our primary interest is in the heat kernel asymptotics. They govern 
the short time asymptotics of quantum mechanical propagators, the ultra violet 
divergences and the large mass expansion in quantum field theory. 

The heat kernel coefficients for singular potentials cannot be obtained in gen- 
eral as limiting cases of smooth configurations. This is clear already from the 
fact that the heat kernel expansion for a smooth potential contains powers of the 
potential taken at the same point. Such expressions become ill-defined in the 
delta-function limit. In our previous paper [TBI] we have observed a surprising 



property that sometimes the ultra violet behavior of a system with 5-function 
potential is less singular than that of corresponding "smooth" system. 

Let us proceed with basic definitions. Let M be a smooth Riemannian man- 
ifold of dimension m. Let E be a smooth closed submanifold of co-dimension 1. 
Let V be a vector bundle over M. Let E and V be endomorphisms of V and V|s 
respectively. In more "physical" notations, E and V are matrix valued functions 
bearing spin and internal indices. In the present paper we study the heat kernel 
expansion for the operator 

D = -(V 2 + E{x) + 6zV(x)) = D - 6xV{x) . (1) 

Let dx and dy be the Riemannian volume elements on M and E. We normalize 
the function 5^ in such a way that for any smooth function / 

/ 5 s f(x)dx = [ f(y)dy. (2) 

We adopt the following short hand notations for the integrals 

/ dxF(x) = {F}[M], [ dyF(y) = {F}[Y]. (3) 

We can choose the coordinates in such a way that in the vicinity of E the 
metric has the form 

g ij dx i dx j = (dx m ) 2 + g ab dx a dx b . (4) 
The second fundamental form of E is L ab = \d m g a i,. We suppose that x m = on 



E. Our notations are the same as in Ref . [J14jJ . Rijki are the components of the 
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Riemann curvature tensor. With our sign conventions, Run is negative on the 
standard sphere in Euclidean space. The Ricci tensor p and the scalar curvature 
r are given by 

Pij ■ Rikkj T Pii Rikki- (5) 

Let p 2 := Pij Pij and R 2 := RijkiRijki be the norm of the Ricci and full curva- 
ture tensors. Let fly be the endomorphism valued components of the curvature 
of the connection on V. In physical language, Qij is the field strength for the 
Yang-Mills and spin-connections. Let ';' denote multiple covariant differentiation 
with respect to the Levi-Civita connection of M and let ':' denote multiple tan- 
gential covariant differentiation on the boundary with respect to the Levi-Civita 
connection of E; the difference between these two is measured by the second 
fundamental form. For example, f ;jj = f. mm + f. aa = f. mm + f :aa - L aa f. m . 

A mathematically rigorous way to define the spectral problem for the operator 
D is to replace it by the spectral problem for D for x E supplemented by the 
conditions on E (see e.g. fl|]) 

0(-O) = 0(+O) (6) 
V m 0(-O) - V m 0(+O) = V<P(0) (7) 

with the short hand notation 0(±O) = lim x .m^_|_ <j)(x). 

The heat kernel K(x, y; t) is a solution of the heat equation 

(d t + D )K(x,y;t) = (8) 

with the initial condition 

K(x,y;0)=5(x,y), (9) 

which satisfies the matching conditions on E following from (Q) and (0). 
We are interested in the integrated heat kernel 

K(f, D; t) = {f(x)K(x, x- t)}[M] = Tr(f exp(-tD)) . (10) 

On manifolds with boundary with local boundary conditions there is an asymp- 
totic expansion as t — > +0 

K(f,D;t) = Y,^(f,D)&- m ^ 2 , (11) 

n=0 

where the coefficients a n (f, D) are volume and surface integrals of local invariants. 
The existence of the asymptotic expansion (|TTD is usually considered as granted. 
There are however important exceptions (besides "genuine" pseudo differential 
operators, the square root of Laplacian for instance), such as the boundary value 
problem for spectral boundary conditions and 5-function potentials with point- 
like support on manifolds with dimension m > 2. In such cases lnt terms can 
appear in the asymptotic expansion ]T5], |T6|, [T7|]. To the best of our knowledge, 
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the existence of the expansion ( |TT| ) for the problem considered here has never 
been stated before. 

This paper is organized as follows. In the next section we derive an integral 



equation for the heat kernel and show validity of the asymptotic expansion (|TT|) . 
In section 3 we calculate the heat kernel asymptotics for the particular case when 
E is a sphere in R m . In section 5 we derive explicit expressions for the heat kernel 
coefficients a n (f, D), n < 5 for the most general form of the operator D. To this 
end we use the particular case calculations of section 4 and functorial properties 
of the heat kernel. 

2 General structure of the heat kernel 

To study the general structure of the heat kernel expansion we use an integral 



equation similar to that proposed by Gaveau and Schulman [|18j for the one- 
dimensional (^-potential: 

K(x,y;t) = K (x,y;t) + f ds f dz K (x, z;t - s)V{z)K{z,y; s) (12) 

Jo Jt, 

where K (x,y;t) denotes the heat kernel corresponding to the operator D with 
7 = 0. 

The equation ([121) has a solution in the form of the power series in 7: 

°° pt ps n pS2 p p 

K(x, y; t) = K (x, y; t) + y~] / ds n / ds n _i . . . / ds x dz n ... dz 1 

n=l ^° ^° ^° ^ s ^ s 



xK (x, z n ; t - s n )V{z n )K Q {z n , z n _ x \ s n - s n _i) . . . V(z 1 )K (zi,y; si) . (13) 

The equation (0) can be obtained formally as a limiting case of the smooth 
potential. Instead of investigating this limiting procedure we prefer to check di- 
rectly that the heat kernel defined by ([12]) satisfies the heat equation (||) with the 
initial condition (Q) and boundary conditions which follow from (|6|) and (|7|). The 
initial condition (|9]) is evident from the eq. (0). Only the first term contributes 
at t — if x, y ^ E. The heat equation (H) can be checked by a direct calculation. 
The first of the matching conditions (|6]) just expresses the fact that the heat 
kernel Kq(x, y; t) is smooth enough. The second condition (0) is a bit less trivial. 
Let y ^ E. Then the following sequence of transformations holds: 



JY, 



-lim / dx m V 2 m / ds / dzK (x,z;t-s)V(z)K(z,y;s) 



lim J dx m Dl J ds j dzK (x,z;t- s)V(z)K(z,y; s) 
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lim / dx m ds dz (d s K (x, z;t - s))V (z) K (z , y; s) 
^°J-e Jo Jt. 

lim J dx m ds J dz K (x, z;t - s)V(z)D y K( Zl y; s) 



+ lim / dx m / dz[8(x,z)V(x)K(z,y;t) - K (x,z;t)V(z)6(z,y)] 

J-e Jt 

= V(x a ,0)K((x a ,0),y;t) 

The main advantage of the representation (|T3|) is that the small t behavior of 
K(x,y;t) is defined through the small t behavior of Ko(x,y;t) which is known 
in some detail. To simplify the notations we do not write down explicitly the 
volume elements here and the parallel transport matrices later on. We also drop 
all matrix indices. 

We are interested in the integrated heat kernel K(t) = J dxK(x,x;t), where 
we put the smearing function / = 1 for simplicity. The integration over x can be 
performed by using the equation 

dxK (x,y 1 ;T 1 )K (y 2 ,x;T 2 ) = K (y 1 ,y 2 ;T 1 + r 2 ) , (14) 

M 

which follows from the evident operator identity e~ TlD °e~ T2D ° = e - (" rL+T2 )' Do . We 
have 

K(t) = K (t) + / ds n / ds n ^i . . . dsi / dz n . . . dz x 

xK Q (zi,z n ;t + si - s n )V(z n )K ( )...V{z x ). (15) 



It is instructive to calculate explicitly several first terms of the expansion (|T5|) . 
In the linear order in V one immediately gets 

K 1 (t)=t / dzK {z,z;t)V{z) . (16) 
Jt, 

This equation means that the liner order in V can be obtained from the asymp- 
totic expansion for the heat kernel in which a smooth potential is replaced by the 
singular one. Such simple relation does not hold at higher orders in V. 

In analyzing the order V 2 contributions we suppose that D is just the stan- 
dard scalar Laplacian in R m , neglect derivatives of V and suppose that S is flat. 
However, we turn to a bit more general case allowing E to be of dimension m — k. 
In this particular case Ko(x,y;t) = (47rf)~ m//2 exp(— (x — y) 2 /At) and 

r rt pS2 

K 2 (t) = (4n)-—t— / dzV(z) 2 / ds 2 / dsi((t-s 2 + Si)(s 2 + Sl )) 2 . 

Jt, Jo Jo 

(17) 
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For k > 1 the integrals over s are divergent. This shows that the proposed 
method can not be extended in particular to 8(r) potentials in R m with m > 2 
for which the expansion (pTTj) is known to break down 0). Return to the 
subject of our study, k — 1. All integrals are easily calculated giving 



x 2 m = — * 1 / cm^) 2 + . . . 



18) 



where we have omitted all higher order terms. 

To calculate the other terms in the expansion (|I5 ) one can use the following 



strategy. Let us substitute the small s asymptotic expansion for the K (x,y; s) 

*><*.Vi •) ~ "-^^^(Af^y) + o*S»(*,v) + ■■■). (19) 

where a(x, y) is the half square of the geodesic distance between x and y, (x, y) 
are the heat kernel coefficients for the operator D . Contributions of largely 
separated points to (|15|) are exponentially damped. Therefore, we may expand 



A*2i(zj, Zj-i) and a(zj,Zj_i) in Taylor series in (zj — Zj-i). The potentials V(z. 



are to be expanded around certain point, say z\. At the end, the expression under 



the integrals in ( [Top will become a sum of the monomials 



exp 



4(t -si + s n ) " ' 4(s 2 - si) 



i {N \v,...; Zl ) 



where I^ N ' is a local invariant functional of V, geometric invariants and their 
derivatives calculated at the point Z\. Negative powers of (sj — Sj_i) appear due 
to the expansion of the a(zj — Zj-i) in the exponentials. It is easy to see that 
Mj < Nj/2. 

Integrals over Zi, except for the last one over z%, can be calculated with the 
help of the relation 



with positive real parameters a and (3. Note that eq. fl21~|) is just a particular case 



of fllil) when D is the flat space Laplacian. Integrals of even powers of {x — yx) 
with the same exponential weight are obtained by differentiation of ([H]) with 
respect to a. Odd powers of x can be integrated by using the following obvious 



relation 



' R' 



, . a f3yt + ay%\ ( (y 1 - xf (x - y 2 ) 2 \ . 
dX ^ X J eX H to 4^)=°- (22) 



6 



Before integrating over Sj let us introduce rescaled variables Sj = Sj/t. This 
enables us to extract an overall power of the proper time t. Integrals over Sj 
will give just numerical factors. One can easily see that the strongest possible 
singularity of the integrand has the form ((1 — s n + §i) . . . (§2 — si))~ 5 - This 
singularity is integrable. Hence the integral over Sj always exists. After all the 
integrations have been done one obtains L dziI^ N \V, . . . ; Zi) multiplied by a 
numerical coefficient and a power of \fi. 

Generalization of the procedure proposed in this section for the case of non- 
unit smearing function / is obvious. Therefore, we have demonstrated that for 
the spectral problem considered in this paper the asymptotic expansion (|TTD is 
valid where the coefficients a n (f, D) are integrals over M and E of local invariants. 
Volume terms are the same as in the heat kernel expansion for the operator D . 



3 Penetrable spherical shell 

Before calculating the heat kernel expansion for generic form of the operator D 
consider a particular case of the constant potential V with the support on a 
spherical shell of the radius R: 

V(x)6 s = ~6(r-R). (23) 

K 

Here r is the radial coordinate. Let D in (|IJ) be the standard scalar Laplacian 
in R m . After separating the angular variables the eigenvalue equation takes the 
form 

& (m-l)a + m + m-2)>| =|fc , 



where I = 0,1,2,... is the orbital momentum. After the substitution <f>(r) = 
r ( 2 - m )/ 2 -?/;(y) the equation takes the form 

92 i£ + £W = *V. (25) 



^^'y 2 ty ^^ r y ■y* 2 

where the notation v — I + is introduced. 

Using the standard techniques |fT9 , the zeta function associated with this 
operator can be written in the form 

00 

^=Jtw) k ^' (26) 



where K(t) is the integrated heat kernel. The function r(s)£(s) has simple poles 
in s — I — N (N — 0, |, 2, . . .) whose residua are determined by the expansion of 
the heat kernel for t — *■ 0. It reads 

K(t) ~ Yl ant {n - m)/2 ■ (27) 

n>0 
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It can be easily seen that the coefficients are related to the residua by means of 

«n= R< S_r( s )c(s). (28) 

S ~ 2 

These heat kernel coefficients can be obtained by calculating the zeta function 
starting from the differential equation, Eq. (|25|). The zeta function of the oper- 
ator D can be expressed in the form 

oo 

OO „ r, 

C(«) = — E A / dk k ~ 2 Yk lnfl{lk) ' (29) 

1=0 { 

where fi(k) is the Jost function of the scattering problem corresponding to the 
operator D and 

(2l + m-2) (l + m-3)\ 

is the multiplicity of the orbital eigenvalues. This representation, as it stands, is 
valid for 3fts > — . The Jost function reads [H| 

fi(ik) = 1 + aI v (k)K„(k) . (31) 

In Eq. ( p9]) the contribution resulting from the empty space (it is independent 
from a) is dropped. By this reason there is no contribution corresponding to the 
coefficient a there. 

In fact, we need the residua of the function r(s)£(s). These are delivered 
when inserting the uniform asymptotic expansion for k — > oo and I — > oo of 
the Jost function into Eq. fl29|) . The latter can be obtained simply by inserting 
the known expansions of the Bessel functions into (|5T|) . This expansion can be 
written in the form 

ln/,(zA;) = £; *n,t l v- n . (32) 

n,i 

Here, n = 1, . . . , N m , and N m is the number of the highest heat kernel coefficient 
requested to be calculated. The coefficients X n i are numbers, some first given in 
the appendix. The notation t = yl + (k/is) 2 is used. When inserting (|3"2] ) into 
the r.h.s. of Eq. (p9|), the change of variables k —>■ vk can be made after what 
the expression factorizes. The integral can be calculated easily. It reads 

oo 

d »-^- r(1 -^) +i/2) - 



The sum over I takes the form 

oo 

((d,n) = J2D l v- 2s - n . (34) 

1=0 
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The functions ((d,n) can be expressed in terms of Riemann and Hurwitz zeta 
functions. They are shown in the appendix. 

Now the heat kernel coefficients can be expressed in the form 



a » = - f^E^ n )E^W (N=l 3 -,2, 5 -,...,N m ). (35) 

^ 2 n=l i 

The sum over i runs from i — n to i = n + A( f^^— ] — 1) where [. . .] denotes the 
integer part. In this form they can be calculated immediately using one of the 
standard computer systems for analytical calculations. As a results, we obtain 
that ai = in any dimension. Other coefficients read for m = 3: 

—a a 2 

° 2 = V? ° 3 = ¥ 

—a 3 a 4 
° 4 = 12^ ° 5 = 64 

_ -(a 3 (4 + 21 a 2 )) _ a 4 (l + 2a 2 ) 

° 3 2520 ° 7 1536 ' 1 J 



for m = 4: 



—a a 2 



° 2 = X ° 3 32 

-(a 2 (3 -4 a 2 ) 0r) 

° 4 " ~48~ ° 5 " 1024 (37) 

-(a 3 (-8 + 7a 2 )) - (a 2 (135 + 168 a 2 - 128 a 4 ) y^) 

° 6 ~ 3360 ° 7 ~ 3932 16 ' 
for m = 5: 

—a a 2 

° 2 = 120F ° 3 = 48 

-a 3 a 2 (-2 + a 2 ) 

° 4 = 720F ° 5 = ^ (38) 



(a 3 (-24 + 7 a 2 )) -9a 4 + 2a 



Or = == a? 

5040 v^F 9216 



for m = 6: 



—a a 2 Vtt 



«2 = -77T a 3 



64 256 



a 



3 



a 2 (15 -4a 2 ) y^ 



a 4 = a 5 = (39) 

384 8192 _ ; 

- (a 3 (-20 + 3 a 2 )) a 2 (945 - 1160 a 2 + 128 a 4 ) y^ 

° 6 ~ 11520 ° 7 ~ 3145728 
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and for m = 7: 



—a a 2 

a2 = 12O0F ° 3 = 480 

-a 3 a 2 (-6 + a 2 ) 

a 4 = = a 5 = (40) 

720 v^F 5 3840 v ; 

-(a 3 (-76 + 7a 2 )) a 2 (27- 15a 2 + a 4 ) 

a 6 = ,= a 7 - 



50400 v^F 46080 

4 Calculation of the heat kernel asymptotics 

The information on the general structure of the heat kernel asymptotic obtained 
in Sec. 2 can be summarized in the following Lemma. 

Lemma 1.1) Let N u (f) = / ;m ... m be the v ih normal covariant derivative. There 
exist invariant local formulae a nu (y,D) so that: 

a n (f,D) = {a n (f,D ;x)}[M} + { £ N»(f)a n , v (y, £>)}[£] (41) 

0<u<n-l 

2) If we expand a n) „ with respect to a Weyl basis, the coefficients depend on the 
dimension m only through a normalizing constant. 

3) Consider a transformation which changes the sign before the m-th components 
of all vector and tensor fields and reverses the sign of the extrinsic curvature L a t. 
Under this transformation a n>1/ — > [ — \) v a n ,v 

4) In the linear order of V 

a n (f,D) = { V (z)-^--a n (f,D )}[E}. (42) 

Proof. First assertion is now evident. Volume terms {a n (f, D ; x)}[M] are given 
in the Appendix B. One can also observe that coefficients before monomials con- 
structed from the geometric invariants depend on the dimension m only through 
the factor (47r)~ m / 2 . A more simple way to prove the second assertion is to con- 
sider product spaces M = S 1 x M± and S = S 1 x S 1 exactly repeating the 
corresponding proof for manifolds with boundaries ||. Assertion 3) follows from 
the fact that we can repeat all calculation of section 2 with the replacement 
x m — > —x m . The last statement of Lemma 1 is just a trivial generalization of 
the equation fll~6|) for non-unit smearing function /. Indeed, it is sufficient to 



represent the variation on the r.h.s. of fl4~2|) as 



5E(z 

Assertion 4) is now evident. 



Ko(t) = J dx ds f(x)Ko(z, x;t — s)K (x, z; s) . (43) 

M JO 



10 



Now we can determine several first heat kernel coefficient up to a few yet 
undetermined constants. 

Lemma 2. There exist universal constants c^i, . . . , 05,10 such that 

a (f,D) = a (f,D ) 
ai (f,D) = 

a 2 (f,D) = a 2 (f,D ) + (4rry m / 2 {fV}m 
a 3 (f,D) = (47r)-( m - 1 )/ 2 {c 3 , 1 /\/ 2 }[S] 

a 4 (f,D) = a 4 (f,D ) + (An)- m ^{c 4A fV 3 + ^frV + fEV 
11 1 

+ ^fV:aa ~ —f:mVL aa + -f ]mm V}[T] 

00 o 

a 5 (f,D) = (^-^-^{c^fV^ + c^frV' + c^p^V 2 (44) 
+c 5A fV 2 E + c^fV 2 L 

aaLbb + c 5fi f'V 2 L ab L ab + c 5t7 fV :aa V 
+c 5 , 8 fV; a V a + c 5<9 f ]m V 2 L aa + c 5<w f ]mm V 2 }[T] 

Proof. According to Lemma 1.1) and 2) any coefficient a n contains all local 
invariants of appropriate dimension. Some of the invariants, like e.g. f'VL aa , 
f- m V etc are ruled out by Lemma 1.3). All terms linear in V are determined by 
Lemma 1.4). 

Lemma 3. c 3j i = |, c 4j i = |, c 5i i = ^, c 5i4 = |, c 5i5 = — ^g, c 5i6 = j^. 
Proof. The coefficients C3 ; i, c 4j i, C5 ; i, c^^ and 05^ are easily calculated using the 
example of the delta-potential on the sphere of the previous section. To calculate 
the coefficient C5 4 consider the case then E = el is a constant proportional to 
the unit matrix. In this case K(f;t) = K(f;t)\ E=0 exp(te). This immediately 
gives c 5)4 = c 3> i = |. 

Several more universal constants can be calculated by a reduction to Dirichlet 
and Neumann boundary value problems. All necessary definitions and explicit 
expressions for the heat kernel coefficient for that problems can be found in the 
Appendix B. 

Lemma 4. 1) Let M = E x [—a, a]. Let V m = d m and let all geometric 
invariants and the smearing function f be symmetric under x m — > — x m . We 
suppose that f and a sufficient number of its derivatives vanish at x m = ±a. 
Then a n (f, D) = a n (f, D , B~) + a n (f, D Q , B + ) where the heat kernel coefficients 
on the right hand side are calculated on E x [0, a], and S — \V . 
Proof. Since reflection of the mth coordinate commutes with D one can subdivide 
the spectral resolution in the two sets (A^, 4>±) with normalized eigenfunctions 
4>±{— x m ) = ±4>±(x m ). Then the heat kernel becomes 

K(f;t) = K.(f;t) + K + (f,t) (45) 

where 

K±(M= ! dxf(x)J2eM-^±)M*) 2 
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f dx m [ dzf(x m ,z) Vexp(-tA ± )( v / 20± 
Jo Jn N 



x m ,z)) 2 . (46) 



Now we observe that v2(fi± are normalized eigenfunctions of the operator Do on 
X x [0, a] satisfying Dirichlet and Neumann boundary conditions 

<P_\ xm=0 = {V m +-V)<f>+\ x m =o = 0. (47) 

Assertion of Lemma 4 follows immediately. 

Under the conditions of Lemma 4 L a b, p mm and f- m vanish identically We 
can calculate only the coefficients of frV 2 , fV :aa V, fV :a V :a and f ;mm V 2 . 
Corollary c 5 , 2 = ^, c 5J = ±, c 5j8 = ^ c 5fW = ^ 

The rest of the universal constants can be calculated by using the conformal 
properties of the heat kernel which are exactly the same as for the usual boundary 
value problem |20|| . 

Lemma 5. If D(e) = e~ 2 f € D, then j-\ e=Q a n (l, D) = (m — n)a n (f, D). 

Under the conformal transformation the metric g acquires a multiplier e 2 ^ . 
V is transformed to e~^ e V. Basic geometric quantities transform as 

(— \e=0^)ij = 5ikf;j + 5jkf;i — Sijf-k 

( — \ e =oL)ab = —$abf;m — fL a b 

j e \e=o(E) = -2fE + ~(rn-2)f. u 

^-| e =o(r) = -2/r + 2(l-m)/.« 
de 

-r\e=o(Pmm) = ~2j ' Pmm ~ f;aa + (1 ~ m)f;mm , (48) 

de 

where T is the Christoffel connection. We need the following conformal relations 
j e \e=oV 2 r = -4fV 2 r + 2(1 - m )[f. mm V 2 + 2f(V aa V + V a V a ) 

L:aaf;mV ] 

" oV 2 P mm = -4V 2 Pmm f-2f(V aa V + V a V a ) + L aa Vf ym 



de 



d , 

e=0 



de 



+(1 - m)V f. mm 
V 2 E = -AV 2 Ef + (m-2)f(V :aa V + V.. a V a ) 

- X -{m - 2)L aa f. m V 2 + l(m - 2)/ ;mm U 2 



f e U=oV 2 L 2 aa = -W 2 L 2 aa f-2(m-l)V 2 L aa f, 
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— \ e =oV 2 L a bL a i ) — —AV 2 L a bL a bf — 2V 2 L aa f. m 
de 

^- e U=oV aa V = -AV aa Vf-(m-3)(V aa Vf + V a VJ) 

^-\e=oV a V a = -4V :a VJ + 2(V aa Vf + V a VJ). (49) 
de 

To obtain the relation (|49| ) we used integration by parts. Let n = 5. By collecting 
the terms with V aa Vf and V 2 L Ma f. m we obtain 

= 4(1 - m)c 5i2 - 2c 5i3 + (m - 2)c 5i4 - (m - 3)c 5i7 + 2c 5j8 

= -2(1 - m)c 5i2 + c 5j3 - -(m - 2)c 5)4 - 2(m - l)c 5)5 - 2c 5j6 - (m - 5)c 5j9 . 

Solving these equation one obtains 05,3 = 1/192 and 05,9 = —5/384. Below we 
collect for convenience of the reader all universal constants of Lemma 2 

c 3 ,i - -, c 4jl - -, c 5il - -, c 5i2 - -, C 5 ,3 - — , c 5 , 4 - -, (50) 

1 11 5 5 1 

° 5 ' 5 ~ ~256' C5 ' 6 ~ 128' ° 5 ' 7 ~ 24' Cs ' 8 ~ 192' ° 5 ' 9 ~ ~384' Cs ' 10 ~ 64 ' 

5 Conclusions 

In the present paper we have studied the heat kernel expansion for for an ope- 
rator of Laplace type in the presence of semitransparent boundaries. We have 
determined the general form of the asymptotic expansion. Namely, we proved 
the validity of the asymptotic series (fill) . We have calculated explicitly several 
first terms of the expansion for the most general operator of Laplace type and 
arbitrary boundary potential. We believe that this is the most complete study 
done in this field so far. 

Our methods of deriving the heat kernel coefficients admit extensive cross- 
checking. Most of the universal constants can be calculated by at least two 
independent methods. If needed, one can calculate the higher coefficients as well. 
As possible generalizations we can suggest the 5' potentials or even general four- 
parameter family on matching conditions |IJ on the hypersurface S. Another 
possible development of the present results can be the renormalization of quan- 
tum field theory in the presence of singular interactions [21]. We believe that 



semitransparent boundaries provide a more adequate framework for the renor- 
malization than "abrupt" boundary conditions of Dirichlet or Neumann type. 



For most recent work on renormalization with singular potentials see W7 . 
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Appendix A 

Some first non zero coefficients X n ^ appearing in (|32]) read: 



a —a 2 a a 3 -3 a 5 a 

— , Aoo = , X = -|- , As 5 — , Ai 7 = , 

2 ' 8 ' 16 24 ' 8 ' 16 



-5 a 2 
32 ' 

-. . 1085a 5a 3 

^■5,5 — T^F + ~ + T?7T' "^5.7 — — A X 5i g — — h -q^-; 



Xaa ~ 32 


a * X 

~ 64' Xa > 


3 a 2 
6 ~ 16 


, A 4j8 — 


27 a a 3 a 5 
" 256 ' 64 ' 160' 


^5,7 — 


145 a 
64 


3a 3 
32 ' X 


-^5,11 


-693 a 


; -^5,13 " 


1155a 


64 


256 


functions defined in Eq. (R3) 


are 





C(2,n) 


= 2C R (2s + n) 




C(3,n) 


= 2C H (2s + n-l 




C(4,n) 


= ( K {2s + n-2) 




C(5,n) 


= |(C H (2 S + n- 


1 

3i 2 


C(6,n) 


= |(CR(2s + n- 


4)- 


C(7,n) 


= |(C H (2 S + n- 


«5 



lc H (2s + n-l;i)) 
(2s + n -2)) 

5 CH(2s + n _ 3 .l ) + ^ CH(2s + n _!.!)) 
where Cr and Ch are the Riemann and Hurwitz zeta functions correspondingly. 



Appendix B 

In this appendix we give expressions for the heat kernel coefficients for the Dirich- 
let and Neumann boundary value problems. Let M be a smooth compact Rie- 
mannian manifold with smooth boundary DM. Let S be an endomorphism on 
VdM and let ;m be a covariant derivative of with respect to inward unit nor- 
mal. We define the modified Neumann boundary operator B + and the Dirichlet 
boundary operator B~ by 

B + (f> ■= + S<j>)\ dM B-</> := <f>\ 9M (51) 

We set S = for the Dirichlet boundary conditions to have a uniform notations. 
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We need only the case of totally geodesic boundary (L ab = 0). We drop 
certain boundary invariants which vanish under conditions of Lemma 4. Several 
first heat kernel coefficients are p0|, [III p2l 



aoC/.A^) = (47r)-™/ 2 Tr(/)[M] 
OiC/.D,^) = ±^(47r)-( m - 1 )/ 2 Tr(/)[9M] 

o a (/,D,B ± ) = (4vr)- m/2 -Tr{(6FE + Fr)[M] + 12/5[aM]} 

6 

a 3 (/, = ±^(47r)-( m - 1 )/ 2 Tr{/(96E + 16r + 1925 2 ) + 24F ;mm }[9M] 



a 4 (/, -D, S ± ) = (47r)- m/2 — Tr(/(60£ ;fcfc + 60r£ + 180£ 2 + 30fi 2 + 12r ;fcfc 



360 

+5r 2 - 2p 2 + 2i? 2 ) [M] + {f{720SE + 120Sr + A80S 3 (52) 
+1205 ;aa ) + 120 f ; mmS)[dM]) 

a 5 (/, = ±^(4vr)-( m - 1 )/ 2 Tr{/(360E ;mm + 1440E ;m S + 720E 2 

+2880ES 2 + 1440S 4 + 240£ :aa + 240r£ + I20tt ab tt ab 
+20r 2 - 8p 2 + 8R 2 + 480rS 2 + 960S.. aa S + 600S :a ^ :a ) 
+/ ;mm (360£ + 360S 2 + 60r) + A5f, mmmm }[dM] . 

On manifold without a boundary one should keep volume contributions only 
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